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Abstract 

The evolution of quadrupole and octupole shapes in Th isotopes is studied in a fully microscopic 
framework based on nuclear Density Functional Theory. The constrained potential energy maps 
and observables calculated with microscopic collective Hamiltonians, indicate the occurrence of a 



j^ , simultaneous quantum shape phase transition between spherical and quadrupole-deformed prolate 

shapes, and between non-octupole and octupole-deformed shapes, as functions of the nucleon num- 
ber. ^^^Th is predicted closest to the critical point of the double phase transition. A microscopic 
mechanism of this phenomenon is discussed in terms of the evolution of single-nucleon orbitals with 
deformation. 
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Quantum phase transitions (QPTs) present a very active field of researcli in condensed 
matter, as well as in mesoscopic systems as, for instance, atomic nuclei, molecules, and 
atomic clusters. Nuclear QPTs, in particular, correspond to shape transitions between com- 
peting ground-state phases induced by variation of a non-thermal control parameter (number 
of nucleons) at zero temperature l|-|5|] . Most experimental and theoretical studies of first- 
and second-order nuclear QPT have considered quadrupole collective degrees of freedom, 
either for axially- symmetric deformed shapes |6l-ll0| or triaxial shapes ll|, |l2| . Much less 
studied, although potentially very interesting, are transitions related to stable or dynam- 



ical octupole shapes (reflection-asymmetric, pear-like shapes). They have been ana 



employing phenomenological geometric models of nuclear shapes and potentials 
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yzed 



algebraic mode 
functionals 19 



s of nuclear structure 



17|, 



and, more recently, microscopic energy density 



20| . In this Letter we report the first fully microscopic study of a simultane- 



ous double QPT between spherical and quadrupole-deformed prolate shapes, and between 
non-octupole and octupole-deformed shapes in Th isotopes, located at ^^^Th. 

Studies of shape transitions usually start with the calculation of potential energy sur- 
faces (PESs) as functions of collective deformation variables. The microscopic analysis of 
quadrupole and octupole shapes presented in this Letter is based on nuclear covariant density 
functional theory (CDFT) 2l|-|23|. which has successfully been applied to the description of 
a variety of structure phenomena all over the chart of nuclides. 

In the present analysis the relativistic mean-field (RMF) implementation of the CDFT 
framework is employed using the functional PC-PKl 2^, and with pairing correlations 
treated in the BCS approximation. The microscopic PESs are obtained by performing con- 
strained RMF+BCS calculations, with constraints on both quadrupole and octupole mass 



moments 



25 



26|. To study the occurrence of possible phase transitions, one must analyze 



the behavior of observables that can be related to order parameters (equilibrium deforma- 
tions, ground-state charge radii, excitation energies, electromagnetic transition rates, etc.). 
Observables that characterize low-lying collective excitations associated with quadrupole and 
octupole degrees of freedom are determined by the eigenstates of the corresponding gener- 
alized collective Hamiltonians (CH), with deformations as dynamical collective coordinates. 
The self-consistent solutions of constrained RMF+BCS calculations: single-particle wave 
functions, occupation probabilities, and quasiparticle energies that correspond to each point 
on the energy surface, are used to calculate, in the cranking approximation, the parameters 



that determine the dynamics of the CH: the collective potential and inertia parameters as 
functions of collective coordinates. The diagonalization of the Hamiltonian yields the^ exci- 
tation energies and collective wave functions that are used to calculate observables 
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FIG. 1: (Color online) Self-consistent RMF-I-BCS triaxial quadrupole energy surfaces in the /32-7 
plane (0 < 7 < 60*^) (upper panels), and axially-symmetric quadrupole-octupole energy surfaces 
in the /32-/53 plane (middle panels), for 220,224,228rji]^^ rpj^g contours join points on the surface with 
the same energy, and the separation between neighboring contours is 0.5 MeV. Probability density 
distributions in the 132- f^s plane for the ground states 0^ of 220,224,228rp]^ (lower panels). 

The self-consistent energy surfaces of 220,224,228'j]-^^ calculated in the constrained 
RMF-I-BCS model with the PC-PKl density functional in the particle-hole channel and 
a separable pairing force 29J in the particle-particle channel, are plotted in the panels of the 
upper two rows of Fig. [1] The contours join points on the surface with the same energy; suc- 
cessive contours differ in energy by 0.5 MeV. The upper panels display triaxial quadrupole 



PESs in the /32-7 plane. One notices the rather rapid transition from the spherical ^^°Th to 
the well-deformed prolate shape of ^^^Th. An interesting feature of the isotopic evolution 
is the flat prolate minimum in ^^^Th, that extends in the interval < /3 < 0.2 of the axial 
deformation parameter. Flat-bottom potentials allow fluctuations in the deformation pa- 
rameters and, therefore, indicate a possible phase transition. The quadrupole PES of ^^^Th 



is reminiscent of those that we analyzed in the A^ f^ 90 region [30l432j] , in connection with 
first-order shape phase transitions between spherical and prolate axial shapes, first studied 
by lachello using the square- well X(5) model |8|. 

If one considers ^^^Th as a system at the X(5) critical point of a first-order shape phase 
transition then, noticing that the PES is rather steep in 7 and in accordance with the X(5) 
model, the dependence of the PES on the two deformation parameters can be decoupled and, 
instead, the axially-symmetric PES in the /32-/33 plane analyzed. This is shown in the middle 
row of Fig. [H where we plot the PESs of 220,224,228rpj^ calculated using the RMF-^BCS model 
with constraints on the expectation values of the quadrupole moment {Q20), and octupole 
moment (Qso)- The region of actinides with A^ ^ 134 is characterized by pronounced oc- 



tupole correlations [33[ and, indeed, our calculation shows a rapid shape transition between 
A^ = 130 to A^ = 138, from non-octupole to pronounced octupole deformations. The PES of 
the transitional nucleus ^^^Th is rather soft with respect to the octupole deformation in the 
region P2 ~ 0.16. This is consistent with the fact that in the Th isotopic chain the lowest 
excitation energy of the state 3~ is found in ^^'^Th 3^]. The rapid transition and critical 



behavior can also be clearly seen in the probability density distributions in the /32-/33 for the 
ground states 0^ of 220,224,228rj]^^ plotted in the lower panels of Fig. [H The 0^ states are 
eigenstates of the axially-symmetric quadrupole-octupole vibrational collective Hamiltonian 
(2DCH), with parameters determined by the self-consistent constrained RMF-I-BCS calcula- 
tions. The probability densities of ^^°Th and ^^^Th isotopes peak at the spherical and stable 
quadrupole-octupole deformed shapes, respectively, whereas for ^^^Th it is rather extended 
along both (32 and /^a. 

To disclose quantitative signatures of possible shape QPTs, we investigate the dependence 
of observables that can be related to order parameters as functions of the control parameter 
- nucleon number. A critical point of a QPT is characterized by a sudden change in the order 
parameter, even though one expects that in systems with a finite number of particles the 
transition is, to a certain extent, smoothed out. In Figure [2] we display the evolution with 
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FIG. 2: (Color online) Evolution of the energy ratios E[^X)IE[2X) (a), B{E2-2f -^ Of) values 
(b), excitation energies of 0^ states (c), expectation values of the octupole deformation parameter 
(fSs) in the ground state Of (d), B{E3;K'^ = Q- ^ K'^ = 0+) (e), and the isotope shifts of the 
ground-state charge radii: {r'^)gs{A + 2) — (r^)gs(A) (f), as functions of the neutron number in 
Th isotopes. Microscopic values calculated with the 5D and 2D collective Hamiltonians based on 
the PC-PKl density functional are compared to available data. The inset shows the comparison 
between the expectation values of the quadrupole deformation parameter (/32) calculated with the 
5DCH and 2DCH. 



neutron number of several quantities that are directly computed using the excitation ener- 



gies and collective wave functions obtained with the 5DCH (for three-dimensional rotations 
and P2-I quadrupole vibrations) and the 2DCH (for axially-symmetric quadrupole-octupole 
vibrations). The parameters of the two Hamiltonians are completely determined by the 
self-consistent PESs calculated with the functional PC-PKl (cf. Fig. [T]). Figure [2] displays 
the isotopic dependence of the energy ratios E{4:f)/E{2f), B{E2; 2f — t- 0^^) values, excita- 
tion energies of 0^ states, the expected octupole deformations (/^s) of the ground state 0^, 
B[E3] K'^ = 0^—7' K^ = of) values, and the isotope shifts of the ground state charge radii: 
(r^)gs(A + 2) — (r^)gs(A), for a sequence of Th nuclei. The calculated values are compared to 
available data. In general, the agreement between the theoretical results and the data is very 
good, especially considering that the calculation does not involve any additional parameter 
beside those used in the calculation of the self-consistent PESs. 

An important result shown in Fig. [2] is that all the considered quantities present pro- 
nounced discontinuities at ^^^Th. The panels on the left point to a phase transition between 
spherical and quadrupole-deformed prolate shapes, whereas the right panels reveal a phase 
transition from non-octupole to octupole-deformed shapes. For both shape QPTs ^^^Th 
appears to be closest to the critical point. An exception is the energy ratio E{A'f)/E{2t) in 
^^^Th, which experimentally is close to the value of 2.9 predicted by the X(5) model J8| of 
the quadrupole QPT, whereas a much smaller value is calculated by the quadrupole 5DCH. 
However, as shown in the inset, this can simply be due to missing octupole correlations in 
the 5DCH calculation. The inset compares the expectation values of the quadrupole defor- 
mation parameter {P2), calculated with the 5DCH and 2DCH. We notice that, in particular 
for ^^"^Th, the quadrupole-octupole Hamiltonian predicts a considerably larger value of (/32), 
which would then translate in a larger value of the ratio E{Af)/E(2f). 

The quantities determined by the evolution of quadrupole-octupole correlations (panels 
on the right of Fig. [2]), display a sudden change at ^^^Th and then appear to saturate for 
heavier Th isotopes. We have, in fact, performed calculations up to A = 234 and found that 
the PESs in heavier nuclei become extended and soft with respect to P^, but with almost 
constant P2 ~ 0.25 — 0.3. Consequently, the excitation energy of the K'^ = 0]~ vibrational 
state predicted by the 2DCH increases with mass number. This is consistent with the 
increase of the excitation energy of the experimental lj~ state (bandhead of K'" = 0~ band), 
which has been interpreted as a signature of the transition from stable octupole deformations 



to octupole vibrations |l3l-ll6|. 
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A microscopic picture of the softness of the potential with respect to both /32 and P^, 
and the related phenomenon of QPTs in ^^^Th, emerges from the dependence of the single- 
nucleon levels on the two deformation parameters. In the upper and middle panels of Fig. [3] 
we plot the single neutron and proton levels of ^^^Th along a path in the P2-P3 plane. The 
path follows the quadrupole deformation parameter P2 up to the position of the equilibrium 
minimum P2 = 0.16, with the octupole deformation parameter kept constant at zero value. 
Then, for the constant value P2 = 0.16, the path continues from /^a = to /^s = 0.3. In the 
lower panel of Fig. [3] we also show the evolution of the neutron and proton pairing energies, 
which reflect the level density around the Fermi surface (dash-dotted curves). 

The extended quadrupole minimum in ^^^Th (upper panel of Fig. [T]) can be attributed to 
the wide region of low proton- level density around the Fermi surface, similar to the density 



of neutron levels in the rare-earth region at A^ ^ 90 for the X(5) QPT 30N32|. In the 
present case the Fermi surface for Z ~ 90 is far above the Z = 82 shell closure, but for small 
quadrupole deformations the splitting of the intruder state is not large enough to enhance 
the density of states. On the other hand, the single-neutron level density is determined by 
the A^ = 126 closure in lighter Th isotopes, and by the splitting of the intruder state J15/2 
in the heavier isotopes. Therefore, the lighter and heavier isotopes display nearly spherical 
and well-deformed prolate shapes, respectively. The emergence of soft octupole deformation 
(middle panel of Fig. [1]) is determined by the evolution of levels along the /^s-path, that is, 
by the coupling between states of opposite parity that originate from the spherical levels 
gQ/2 and J15/2 for neutrons, and fj/2 and 213/2 for protons J33|. 

In conclusion, using a fully microscopic theoretical framework based on nuclear CDFT, 
we have analyzed the evolution of quadrupole and octupole shapes in Th isotopes. 
The RMF+BCS calculations of constrained quadrupole triaxial and axially-symmetric 
quadrupole-octupole energy maps have been performed. The constrained RMF-I-BCS solu- 
tions determine the parameters of a quadrupole 5DCH, and a quadrupole-octupole 2DCH, 
that are used to calculate observables related to order parameters. Both the calculated 
PESs and the predicted observables (excitation energies, isotope shifts of charge radii, elec- 
tromagnetic transition rates) indicate the occurrence of a simultaneous phase transition 
between spherical and quadrupole-deformed prolate shapes, and between non-octupole and 
octupole-deformed shapes, with ^^^Th being closest to the critical point of the double QPT. 
A microscopic interpretation of the QPT has been given in terms of the evolution of single- 
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FIG. 3: (Color online) Single-neutron levels (upper panel), single-proton levels (middle panel), 
and neutron and proton pairing energies (lower panel) of ^^^Th, as functions of the deformation 
parameters. Each plot follows the quadrupole deformation parameter (32 up to the position of the 
equilibrium minimum (32 = 0.16, with the constant octupole deformation parameter (3^ = {) (left 
panels). For the constant value (32 = 0.16, the panels on the right display the dependence of the 
single-nucleon energies and pairing energies on the octupole deformation, from (3-^ = Q to (3^, = 0.3. 
The thick dash-dotted (black) curves denote the Fermi levels. 

nucleoli orbitals with quadrupole and octupole deformation parameters. The results are 
in good agreement with experiment and with previous studies of shape transitions in this 
region. However, this is the first study in which both mean-field PESs and order parameters 
have been calculated fully microscopically, without any phenomenological control parame- 



8 



ter. Double shape QPT and related critical-point phenomena could occur in other finite 
quantum systems characterized by quadrupole and octupole collective degrees of freedom. 
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